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ABSTRACT 

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) 
mapped at 11 pc resolution in the CO(l-O) line. Targets were chosen based on a limiting CO flux and peak 
brightness as measured by the NANTEN survey. The observations were conducted with the ATNF Mopra Tele- 
scope as part of the Magellanic Mopra Assessment (MAGMA). We identify clouds as regions of connected CO 
emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decom- 
position parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL oc L~ 2 , 
suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, 
while apparent for the largest emission structures, breaks down when those structures are decomposed into 
smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of com- 
paring two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth 
relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide 
range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar 
object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the 
LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more lumi- 
nous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between 
giant molecular clouds and massive star formation. 

Subject headings: Magellanic Clouds — ISM: molecules — galaxies: ISM — Stars: formation 



1. INTRODUCTION 

Recent surveys of CO in Local Group galaxies have 
achieved sufficient resolu tion to resolve individual molecular 
cloud s (see reviews bv lBlitz et al.ll2007l : iFukui & Kawamura! 
2010). These studies represent an important step towards con- 
necting our understanding of star formation in the Milky Way 
to more distant and less evolved systems. Among the ques- 
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tions that these studies seek to address are whether molecular 
cloud characteristics, and the conditions for their formation, 
follow universal patterns, and what the resulting implications 
for star formation may be. 

The Large Magellanic Cloud (LMC) is the nearest star- 
forming galaxy to our own (d « 50 kpc), yet the distribution of 
its CO emission has been revealed in detail only recently, due 
to the LMC's large angular size and the paucity of millimeter- 
wave facilities in the southern hemisphere. The first complete 
C O map, at a ( smoot hed) resolution of 12', was published 
bv lCohen et al.l (1988), based on observations with the 1.2-m 
Columbia Millimeter- Wave Telescope. Similar spatial cover- 
age was later obtained by the 4-m NANTEN telescope, at a 
resolution of 2.6; result s of the init ial NANTEN s urvey were 
reported by IFukui etail (119991) and lMizuno et all (1200 ll) . and 
results of a se cond, more s ensitiv e survey have been recently 
published by IFukui et all d2008l hereafter Fu08). In addi- 
tion to the large-scale surveys, detailed mapping of individ- 
ual cloud complexes (at ~1' resolution, corresponding to 15 
pc) has been pursued with the Swedish-ESO Submillimetre 
Telescope (SEST) to explore the sizes, dynamical state, and 
molecular and isotopic abundances of individual clouds (e.g. , 
[1991 iKutner eTaT] \l99% Uohansson et all Wm. 
l2003l) . 
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The new data set presented in this paper arises from a natu- 
ral extension of the NANTEN surveys: the targeted mapping 
of known CO complexes in the LMC at improved angular res- 
olution (45", or 1 1 pc at a distance of 50 kpc), adequate to re- 
solve the largest giant molecular clouds (GMCs). These maps 
are a principal outcome of the Magellanic Mopra Assessment 
(MAGMA), a long-term project of CO mapping of both the 
LMC and the Small Magellanic Cloud (SMC) with the Mo- 
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pra Telesc ope. Initial results from MAGMA have bee n pre- 
sented bv lPinedaetail ( 120091) a ndlHughes et alj (120101 here- 
after HulO) for the LMC and iMuller et aLllioid) for the 
SMC. In the current work we present the completed LMC data 
set and three nested catalogs of GMCs, obtained using differ- 
ent assumptions for identifying discrete clouds. We compare 
the cloud properties and scaling relations derived from these 
catalogs with previous work and discuss the evidence for and 
against clouds being in a state of virial equilibrium. Although 
our final catalogs differ from the one used in HulO, the distri- 
butions of cloud properties do not change dramatically, and in 
the present paper we do not revisit the relationships between 
GMC properties and environmental conditions (such as stel- 
lar surface densities, interstellar pressure, and far-ultraviolet 
radiation field) examined in HulO. We expect substantially 
improved ancillary (non-CO) data to become available in the 
near future, at which point we intend to re-examine those re- 
lationships. 

As one of the initial applications of the MAGMA data 
set, we examine the correspondence between CO clouds and 
known or suspected young stellar objects (YSOs). While stars 
are believed to be formed in molecular clouds, the relation- 
ship between molecular gas and star formation has been dif- 
ficult to analyze in external galaxies because of observational 
biases. For example, while Ha is commonly used to trace 
recent star formation, it is susceptible to extinction and there- 
fore the correlation between CO and Ha may break down on 
scales approaching that of molecular cloud complexes. At in- 
frared wavelengths, even deeply embedded young stars are 
detectable if they are sufficiently massive, leading to two re- 
cent attempts to obtain a census of massi ve YSOs in the LMC 
based on Spitzer Spa ce Telescope data (Whi tney et al.ll2008t 
iGruendl &^ Chu 2009). As the reliability and completeness of 
such catalogs improve, it should be possible to measure the 
star formation efficiency and lifetimes of molecular clouds 
based on number c ounts of YSOs rather than integrated fluxe s 
of dust or gas (e.g jlndebetouw et alj2008tlChen et al.l2010t) . 

We emphasize that, while we will generally refer to the 
sources we detect as "molecular clouds," the results presented 
in this paper refer only to the CO-emitting molecular clouds. 
The reliability of CO as a tracer of H2 has been a topic of ac- 
tive discussion (e.g.. lLero v et al. 2007; Bot et al. 2007). Un- 
fortunately, the bulk of the H2 is at temperatures too low 
to produce detectable emission, so it can be probed only 
through absorption or through indirect emission tracers such 
as CO, dust and gamma rays. Far Ultraviolet Spectroscopic 
Explorer (FUSE) observations of H2 in absorption towards 
LMC targets suggest a total diffuse H2 mass of 8 x 10 6 M Q 
dTumlinson et alJ2002l) . < 2% of the LMC's atomic gas mass. 
However, the presence of a significant mass of H2 residing 
in CO-dark envelopes surrounding LMC mol ecular clouds is 
not str ongly excluded by this result, as the iTumlinson et al.l 
(2002) sample contains predominantly sightlines with low ex- 
tinction. Indeed, the ratio of CO emissio n to H2 m ass is 
likely to decrease on larger spatial scales (Israel 2000j), and 
studie s of LMC dust emission at far-infrared wav elengths 
(e.g.. iBernard et af]|2008t iRoman-Duval et al1l2010t) suggest 
the presence of more dust than can be accounted for using 
standard conversion factors for H I and CO emission. A prob- 
lem for these studies is whether to attribute variations in the 
inferred dust-to-gas ratio to variations in dust abundance (or 
emissivity) or to incorrect estimation of gas masses. In any 
case, there is general agreement that CO emission traces re- 



gions where hydrogen is predominantly molecular, although 
it may not identify all such regions. 

We begin with an overview of the survey strategy, the obser- 
vational parameters, and data reduction procedures in Section 
[2] The identification of clouds, their derived properties and 
a comparison with YSO catalogs are presented in Section [3] 
Finally, Section |4]presents a discussion and Section [5] a sum- 
mary of our results. 

2. THE MAGMA SURVEY OF THE LMC 

2. 1 . Survey Design 

The MAGMA survey area was limited in total to roughly 
3.6 deg 2 by the need to staff the observations and complete 
them in a reasonable amount of time. Thus, we chose to map 
regions coincident with, or close to, known CO clouds de- 
tected by NANTEN. The second NANTEN survey (Fu08) de- 
tected 272 clouds, of which 230 were detected in more than 2 
grid positions (spaced 2' apart) and termed "giant molecular 
clouds." Hereafter we refer to these 230 objects as NANTEN 
GMCs, and the larger sample of 272 clouds simply as NAN- 
TEN clouds. When the scope of MAGMA was expanded 
from limited regions in the eastern and central part of the 
LMC to the entire galaxy, we adopted both flux and brightness 
criteria for NANTEN GMCs to be included in our survey. A 
minimum extrapolated CO flux of 1.2 x 10 5 K km s" 1 arcsec 2 
(equivalent to a CO luminosity of 7000 K km s -1 pc 2 at the 
distance of the LMC) was chosen in order to ensure that the 
largest GMCs were included. (Details of the flux extrapola- 
tion used by Fu08 and this paper are given in Section [3~2l ) In 
addition, a minimum peak CO brightness of 1 K km s" 1 was 
chosen to exclude a handful of clouds (Nos. 34, 38, 78, & 
138 in Fu08) which met the flux requirement due to a large 
extrapolated size but were deemed unlikely to be detected by 
Mopra. Our fiducial sample was thus limited to 114 NAN- 
TEN GMCs, although in fact many of the remaining GMCs 
ended up in the MAGMA field because of their proximity to 
survey targets. A detailed compariso n of the MAGMA and 
NANTEN samples is given in Section l2~4l 

2.2. Observations and Data Processing 

Observations were performed with the ATNlFl Mopra 22- 
m telescope from 2005 May to 2010 October. We divided our 
target regions into square 5 '2 x 5. '2 patches whose centers 
were separated by 4 .'75, ensuring significant overlap between 
adjacent patches. Each patch was mapped at least twice in 
the On-The-Fly (OTF) observing mode, in which spectra are 
taken continuously as the telescope is scanned across the sky. 
Scanning is conducted in either right ascension or declination; 
to minimize scanning artifacts we alternated the scan direction 
when repeating a patch. Along each row of an OTF field, indi- 
vidual spectra are recorded every 7", so that the 33" (FWHM) 
telescope beam is oversampled in the scanning direction. The 
spacing between rows is 10", also oversampling the beam. In 
our standard configuration, an OTF map consists of 3 1 rows 
with 44 2-second samples taken per row, so each map con- 
tains about 2700 spectra (counting both polarizations sepa- 
rately). A 20-second off-source spectrum, taken towards a re- 
gion within 1° that shows little H I emission (J 7b dv < 280 K 
km s" 1 , corresponding to Nh <5x 10 20 cm -2 ), is observed be- 
fore each row, and the above-atmosphere system temperature 

16 The Australia Telescope is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO. 
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Figure 1. CO intensity contours from MAGMA, at levels of 2 and 10 K km s -1 , overlaid on an Ha image from MCELS iSmith et al. 1999), shown on a 
logarithmic stretch to bring out fainter emission. The CO image was derived by smoothing the MAGMA cube to 1' resolution and applying a blanking mask 
corresponding to the 3<r contour of a 1.5'-resolution smoothed cube. Multiple small rectangular contours show the regions mapped in CO with Mopra. Large 
rectangular boxes show the regions magnified in Figures |5]and|6] The Ha image has not been continuum subtracted, although bright point sources have been 
removed by a moving boxcar median filter. 



is estimated every 25 minutes using an ambient (hot) load. 
Typical system temperatures for the survey observations are 
between 500 and 600 K; observing is abandoned when system 
temperatures exceed ^1000 K. Fluctuations in system tem- 
perature between hot load measurements are monitored using 
a rapidly switched noise source injected at the front end of the 
receiver. A single OTF map takes about 75 minutes to com- 
plete. After each map, the pointing of the antenna is verified 
and updated by observing the nearby SiO maser R Dor. Prior 
to correction, the pointing errors are typically below 10". The 
final MAGMA CO data cube includes data from 1163 OTF 
maps with approximately 3 100000 individual spectra. Fig- 
ure [U shows the spatial coverage of the MAGMA data set 
in comparison to an Ha image from the Magella nic Clouds 
Emission-Line Survey (MCELS; Smit h et al.Hl999l) . 



Data reduction was performed using the ATNF's MIRIAD, 
Livedata, and Gridzilla software packages. Livedata per- 
forms bandpass calibration using off-source spectra, fits and 
subtracts a linear baseline, and applies additional calibration 
factors based on assumed telescope efficiencies, which we 
derived from an analysis of st andard spectra taken towards 
Orion KL (cf. lLadd et al]l2005h . The efficiencies we have as- 
sumed for CO are listed in Table Q] and represent the "ex - 
tended beam" efficiency rj^ discussed by lLadd et al.l (12005b . 
which includes the effect of coupling to the inner error beam 
for sources larger than ~2'. We determined rj x \, for each stan- 
dard spectrum by dividing the observed peak antenna temper- 
ature by 100 K, which is the correc ted peak CO brig htness 
temperature for Orion KL given bv lLadd etalJ d200l . An 
average value for r] x \, was then calculated for each observ- 
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Table 1 

Adopted Efficiencies for the Mopra Telescope 



Start Date 


End Date 


Vxb 


2005-05-29 


2005-10-10 


0.54 


2005-10-20 


2005-10-29 


0.38 


2006-06-25 


2006-10-27 


0.38 


2007-07-01 


2007-09-30 


0.45 


2008-06-19 


2008-06-29 


0.39 


2008-08-17 


2008-09-30 


0.48 


2009-06-15 


2009-09-30 


0.48 


2010-07-26 


2010-10-10 


0.43 



ing season with a typical rms deviation of — 10%. Note that 
some of the apparent changes in r/ x b likely reflect changes in 
the system temperature scale due to instrument modifications 
rather than changes in the aperture efficiency. Gridzilla takes 
the calibrated spectra and grids them onto a data cube using a 
Gaussian smoothing function of FWHM 33", comparable to 
that of the Mopra primary beam. We weight the spectra by 
the inverse of the system temperature when gridding. Based 
on simulations of the gridding process we estimate that the 
final resolution of our maps is close to 45". 

Changes in the spectrometers available at Mopra between 
2005 and 2006 have resulted in varying spectral resolution 
and bandwidth over the period of our observations. The spec- 
tral resolution varied from 33.6 to 62.4 kHz in 2005, so data 
from that year were resampled by spline interpolation to the 
resolution of 33.68 kHz (0.09 km s" 1 ) available in later years 
using the MIRIAD task MOPFIX before processing with Live- 
data and Gridzilla. Since all of our analysis is performed at a 
much coarser resolution (0.5 km s" 1 ), this resampling has lit- 
tle effect on the final data products. Gridzilla generated a data 
cube with 1599 channels spanning radial velocities of 180 to 
320 km s" 1 in the LSR frame; we subsequently binned the 
data in velocity, reducing the number of channels by a factor 
of 6. To summarize, our final output cube, produced using 
an orthographic projection, has an effective spatial resolution 
of 45", sampled onto a grid spacing of 15", and a channel 
spacing of 0.526 km s" 1 . 

2.3. Integrated Intensity Images 

We explored three different methods to construct CO in- 
tensity images. The first, and most straightforward, is di- 
rect integration over the expected velocity range. This re- 
sults in excessive noise when applied to the entire galaxy, 
although limited application to small regions of the LMC is 
generally more successful [see Figure [2(a)], because smaller 
velocity ranges can be chosen. The second and third meth- 
ods involve integration over a mask that isolates significant 
emission from noise. For the second method, which we call 
the smooth-and-mask method, the original cube is convolved 
with a Gaussian in the spatial domain to achieve a coarser 
resolution (90"), and then a 3a threshold level is chosen in 
the smoothed cube to generate a mask which is applied to 
the original cube. For the third method, which we call the 
dilated mask metho d, w e applied the dilated CPROPS mask 
described in Section [3~2l below. This begins by identifying re- 
gions of high significance (>3cr) and expanding to connected 
regions of lower significance (>2<r). Examples of applying 
the second and third methods to the region around a particular 
GMC (No. 69 in Fu08, also depicted in Fig. 2 of HulO) are 
shown in Figures |2jb) and (c) respectively. As discussed be- 



Table 2 

Comparison of Global CO Fluxes From Various Methods 



Method 


NANTEN data 


NANTEN data 


MAGMA 




(full) 


(MAGMA region) 


data 


Direct integration 


8.2 


6.3 


14 


Smoothed 3a cut 


7.2 


5.8 


4.7 


CPROPS islands 


5.0 


4.7 


3.0 



Note. — Units are 10 7 K km s 1 arcsec 2 . 



low (Section l2~4l . the smooth-and-mask method provides the 
best agreement with the CO flux measured in the NANTEN 
cube. We show the CO intensity derived from this method 
as contours in Figure [TJ On the other hand, the dilated mask 
method, being more conservative, may be preferred in situ- 
ations where the reliability of signal identification is the pri- 
mary concern. 

At the 45" resolution of our final channel maps, the RMS 
brightness temperature in a 0.5 km s" 1 channel ranges from 
0.2 to 0.5 K with a mean value of 0.3 K. For a typical 
equivalent line width of \/2na v — 3 km s" 1 , or 6 spectral 
channels, this corresponds to an RMS integrated intensity 
of 0.4 Kkms" 1 . Thus, for a Galactic Xcp-f actor of 2 x 
10 20 cm- 2 (Kkm s" 1 )" 1 (Strong & Mattox 1996), the ia sen- 
sitivity limit is ~ 1.2 K km s" 1 5 M Q pc~ 2 . This sensitivity 
limit is optimistic as it assumes one can correctly identify the 
mean CO velocity; without such prior information one would 
need to integrate over a wider velocity range, thus increasing 
the map noise. Moreover, if the appropriate value of Xqq is 
larger than the Galactic value, as is commonly assumed, the 
sensitivity limit would scale upward accordingly. 

2.4. Comparison with NANTEN 

The NANTEN CO cube presented by Fu08 contains a total 
flux of 8.2 x 10 7 Kkms" 1 arcsec 2 when integrated directly, 
or 7.2 x 10 7 K km s" 1 arcsec 2 after applying a smoothed (to 
6') 3cr-contour mask0 The corresponding fluxes, when re- 
stricted to the regions surveyed by Mopra, are 6.3 x 10 7 
K km s" 1 arcsec 2 and 5.8 x 10 7 K km s" 1 arcsec 2 respectively. 
Thus, both approaches show that the MAGMA field con- 
tains -80% of the total CO flux observed by NANTEN in 
the LMC. The Mopra CO cube contains a total flux of 1.4 x 
10 8 Kkms" 1 arcsec 2 when integrated directly, or 4.7 x 10 7 
Kkms" 1 arcsec 2 when applying a smoothed (to 1.'5) 3a- 
contour mask. The surprisingly large flux obtained by direct 
integration of the Mopra cube appears to result from system- 
atic errors in the spectral baseline which accumulate when 
summing large numbers of noise channels. We consider that 
measurement unreliable, and conclude that integrating the 
spectra over their full velocity range will not yield reliable 
line fluxes. On the other hand, the fluxes measured by NAN- 
TEN and Mopra using the smooth-and-mask technique are in 
reasonable ag reem ent (to within —20%; see Table[2]i. 

In Section 13.21 we discuss how we used a dilated mask 
to identify regions ("islands") of significant emission. The 
total CO flux summed over these 450 islands is 3.0 x 10 7 
K km s" 1 arcsec 2 . This represents only 64% of the total CO 
flux as determined by the smooth-and-mask technique, sug- 
gesting that a significant fraction of the CO flux lies near or 
below the limit imposed by our CO detection criteria (cf. Fig- 

17 For CO(1-0), 1 K km s" 1 arcsec 2 = 9.6 X 10~ 3 Jy km s" 1 = 3.7 X 10" 23 
WnT 2 . 
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Figure 2. Integrated CO intensity maps for a relatively isolated GMC (No. 69 in Fu08) produced using three different methods: (a) Direct integration of a 10 
km s window centered on the emission; (b) masking using a 3<r threshold applied to a smoothed cube; (c) signal detection with CPROPS as described in the text. 
Contour levels range from 2 to 10 K km s in increments of 2 K km s . The ellipses in each panel (left to right) indicate the extrapolated sizes and orientations 
of structures identified by CPROPS using the "islands", "physical", and "data-based" decomposition parameters respectively; note that these structures have been 
identified in the data cube and not in the CO intensity maps shown. The Mopra beam size (Gaussian FWHM) is represented by the blue circle at the lower left of 
each panel. 
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Figure 3. Left: Map showing as shaded regions the NANTEN CO clouds that do not overlap with MAGMA clouds in position and velocity. These constitute 
41% (112/272) of the NANTEN cloud sample, and primarily lie outside the regions observed with Mopra (shown as blue rectangular contours). Violet contours 
show the NANTEN CO integrated intensity (at levels of 1, 4, 9, and 16 K km s ), with a light gray contour at the periphery outlining the region observed with 
NANTEN. Right: Map showing as red ellipses the locations of MAGMA clouds near the 30 Doradus region in the southeast of the LMC that do not overlap 
with NANTEN clouds in position and velocity. The MAGMA CO intensity map is shown in grayscale at contour levels of 2 and 10 K km s , overlaid with the 
NANTEN CO intensity (violet contours) at levels of 0.75, 3, and 6.75 K km s -1 . 
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ure[2]i. We explored this possibility by relaxing the threshold 
for including emission that is contiguous with strong emission 
(i.e., further dilating the mask). Lowering the threshold from 
2cr to 1 .5ct increases the flux contained with the mask by 26%, 
implying that CO emission is indeed present just below our 
detection limit. For comparison, Fu08 identified significant 
CO emission in the NANTEN cube using a dilated mask as 
well, with a total masked CO flux, restricted to the MAGMA 
observed region, of 4.7 x 10 7 K km s" 1 arcsec 2 . This is 81% 
of the total NANTEN-measured flux over this region, when 
determined using the smooth-and-mask technique, and 75% 
of the NANTEN flux when determined using direct integra- 
tion (Table |2). 

While the NANTEN survey catalogued 272 clouds, we 
stress again that the MAGMA survey covers a much more 
restricted region of sky. For 85 of the "islands" identified in 
the MAGMA cube, there is a one-to-one correspondence with 
a NANTEN cloud. An additional 275 "islands" correspond 
to 72 NANTEN clouds in such a way that more than one "is- 
land" overlaps with a single NANTEN cloud. Finally there 
were a handful of matches (3 islands, 6 NANTEN clouds) 
where two or more NANTEN clouds overlap a single island. 
Removing clouds common to the latter two groups, we find 
a total of 160 NANTEN clouds associated in both position 
and velocity with MAGMA islands, underscoring our earlier 
point that while only 114 NANTEN clouds were originally 
targeted, a larger number were ultimately observed. The ad- 
ditional clouds are mostly lower-luminosity GMCs, although 
4 are considered "small" clouds in Fu08. 

The left panel of Figure [3] shows in color scale the popu- 
lation of NANTEN clouds not observed or not detected by 
MAGMA, while the right panel shows as red ellipses the pop- 
ulation of MAGMA clouds in the "molecular ridge" south of 
30 Doradus which were not detected by NANTEN. The ef- 
fective sensitivity of the two surveys appears to be quite com- 
parable, as evidenced by the fact that most of the NANTEN 
clouds not detected by MAGMA lie outside the MAGMA ob- 
served regions, and most of the MAGMA clouds not detected 
by NANTEN are small clouds at the periphery of NANTEN- 
detected clouds. Of the original 114 targeted NANTEN 
GMCs, only one (No. 46 in Fu08) was not formally detected 
by MAGMA; this can be seen near a2ooo=5 h 4 m , (52ooo=-68°. 
Thus, even though the MAGMA survey is a factor of —2 
worse in brightness sensitivity, this is compensated by the fact 
that most (if not all) GMCs have unresolved structure at the 
NANTEN resolution and therefore appear brighter when ob- 
served with Mopra. 

To summarize, based on the NANTEN data, the observed 
MAGMA field covers a large fraction (—80%) of the total CO 
emission from the LMC, and an even larger fraction (—94%) 
of the CO flux found in catalogued clouds. This is despite the 
much smaller area covered by MAGMA (3.6 deg 2 rather than 
30 deg 2 ). Total CO fluxes over the MAGMA field, as deter- 
mined using the smooth-and-mask method, are consistent be- 
tween the NANTEN and MAGMA cubes. However, the frac- 
tion of the total flux that is found in CO-detected "clouds" is 
lower for MAGMA (64%) than for NANTEN (8 1 %) over this 
same area. This may reflect the poorer brightness sensitivity 
of the Mopra observations, resulting in a significant fraction 
of emission lying just outside the cloud boundaries, as well as 
the fact that the MAGMA fields enclose the NANTEN clouds 
quite tightly, so that the clouds contribute a disproportionate 
fraction of the NANTEN flux in these regions. Indeed, the 
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Figure 4. Distribution of the logarithm of the molecular gas surface density, 
E, (derived from the integrated CO intensity using a constant Xqo) in the 
observed MAGMA field. The black histogram shows values obtained when 
using the smooth-and-mask method to derive the CO intensity, while the red 
histogram shows values obtained using the dilated mask method. A vertical 
dashed line represents our nominal 3cr sensitivity limit of 1.2 K km s 5 
Mq pc~ 2 . The dotted parabola represents a lognormal distribution of surface 
density with a mean value of 16 Mq pc~ 2 and a dispersion above this value 
that is consistent with the data in the black histogram. 

fraction of flux in clouds for the NANTEN survey as a whole 
is only —70%. We therefore stress that the MAGMA survey's 
strength lies in the characterization of the bright CO clouds, 
and that due to limited spatial coverage and sensitivity it pro- 
vides relatively weak constraints on the CO emission flux (let 
alone H2 mass) of the LMC as a whole. 

3. RESULTS 

3.1. Column Density Distribution 

A histogram of CO integrated intensity across the MAGMA 
field can be converted into an H2 column density distribu- 
tion under the assumption of a constant CO-to-FL conver- 
sion factor (Xco)- Figure |4] shows the CO intensity distribu- 
tion, converted to mass surface density £ using Xco = 2.0 x 
10 20 cm" 2 (Kkm s -1 ) , derived from two different methods 
of computing the integrated CO intensity, the smooth-and- 
mask approach and the dilated mask approach. The dilated 
mask values, given by the red hist ogra m, are clearly more 
conservative as discussed in Section [231 and as expected are 
strongly suppressed below the sensitivity limit. A lognor- 
mal function with peak position chosen by eye and a width 
determined by the RMS value above the peak is shown for 
comparison. In the high column density regime, the distribu- 
tion can be reasonably approximated as lognormal, although 
there may be a cutoff at the high end around £ — 200 M Q 
pc" 2 . A lognormal distribution in column density is pre- 
dicted by numeri cal simulations of the turbulen t interstel- 
lar medium (e.g.. iBallesteros-Paredes et alj|201 11) . although 
an excess at low column densities, such as we find in Fig- 
ure |U is o ften seen as well, representing warmer or unper- 
turbe d gas (IWada & Norman 20071 : IBallesteros-Paredes et alJ 
2011). Given the limited sensitivity and biased spatial cov- 
erage of MAGMA, the amount of mass in the low column 
density regime is not well constrained. 

3.2. Cloud Identification 

To identify significant CO emission in the MAGMA 
data cubes, we have used the algorithms presented by 
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Table 3 

Properties of Clouds Derived from CPROPS Analyses 

Islands Physical Data-based 



No. of clouds 


450 


543 


694 


No. of isolated clouds 


243 


222 


190 


Tot. M\ am /MQ (unextrap.) 


7.4 X 10 6 


5.4 x 10 6 


4.1 x 10 6 


Total M lum /M 


1.2 X 10 7 


1.1 x 10 7 


1.1 x 10 7 




0.27 


0.42 


0.55 


Mean(rms) logM| umiM0 


4.00(0.54) 


4.02(0.49) 


3.97(0.44) 


Mean(rms) logi? pc 


1.13(0.27) 


1.09(0.22) 


1.05(0.20) 


Mean(rms) logc^, kms -i 


0.18(0.20) 


0.17(0.19) 


0.14(0.19) 



"Fraction of the total M\ am in identified clouds that is in clouds with masses 
below the completeness limit of 3 X 10 4 Mq. 



iRosolowsky & Leroyl (120061 hereafter RL06), implemented 
in IDL as part of the CPROPS package. CPROPS uses a di- 
lated mask technique to isolate regions of significant emission 
within spectral line cubes, and a modified watershed algo- 
rithm to assign the emission into individual clouds. Regions 
of significant emission within the MAGMA data cubes are 
initially identified by finding volume pixels ("voxels") with 
emission greater than a 3cr threshold across two contiguous 
velocity channels, where the RMS noise a is estimated at each 
sky position using the median absolute deviation of the spec- 
trum. The mask is then expanded to include voxels that have 
Tb > 2a in at least two adjacent velocity channels, and in ad- 
dition can be joined to the initial mask through a path that 
only passes through valid (unmasked) voxels. If a region de- 
fined this way (referred to as an "island") has a projected area 
less than a single telescope beam, or spans a velocity range 
of less than 4 channels, it is assumed to be a noise peak and 
is removed from the mask. Note that HulO used slightly dif- 
ferent masking criteria, requiring a 4a peak and expanding to 
the 1.5 a edge, and adopted a single value of a for the entire 
cube. 

Once "islands" of significant emission have been identified, 
CPROPS decomposes the emission into individual "cloud" 
structures. Moments of the emission along the spatial and 
spectral axes are used to determine the size, linewidth and 
flux of the clouds, and corrections for the finite map sensi- 
tivity and instrumental resolution are applied to the measured 
cloud properties. To explore the effects of tunable parame- 
ters on the results, the decomposition procedure is performed 
under three different sets of assumptions: 

1. "Islands" parameters: CPROPS simply catalogs all 
contiguous emission structures ("islands") in the mask, 
regardless of size (but recall that islands must span 
an area of at least one telescope beam). The size, 
linewidth, and flux of each island is determined by the 
standard RL06 algorithm, involving extrapolation to a 
zero-intensity boundary and correction for finite reso- 
lution of the telescope and spectrometer (see below). 

2. "Physical" parameters: Islands are decomposed into 
clouds, defined as local maxima within a moving box of 
size 5x5 pixels (sal 8 pc) and velocity width 9 chan- 
nels (w4.7 km s" 1 ). These scales are chosen to be sen- 
sitive to structures that observers commonly define as 
GMCs in the Galaxy. The algorithm requires that local 
maxima meet certain distinctness criteria (they must lie 
at least 1 K above the merge level with any other maxi- 
mum, and yield significantly different moments when 



merged with adjacent maxima), and assigns to each 
maximum the surrounding emission that lies above the 
merge level with all other maxima. Sensitivity and res- 
olution corrections are then applied. 

3. "Data-based" parameters: Clouds are identified as lo- 
cal maxima within a moving box of size 5x5 pixels 
and 9 channels, as in the "physical" case. Each local 
maximum is required to lie at least 2a («0.6 K) above 
the merge level with any other maximum. Sensitivity 
and resolution corrections are then applied. 

The results of the three decompositions are summarized 
in Table [3] It is clear that each of these approaches has 
advantages and disadvantages. The "islands" approach ig- 
nores substructure, and defines clouds based on boundaries 
which may depend on the sensitivity of the observations. The 
"data-based" approach identifies all structures which would 
be easily distinguished by eye, but the resulting sizes tend 
to lie near the resolution limit, hindering comparison with 
studies performed at different spatial resolution. The "physi- 
cal" approach is similar to the data-based approach, with the 
main difference being the stricter constrast requirement for 
distinctness (1 K w 3er rather than 2a). We note that the 
decomposition of islands into clouds (whether "physical" or 
"data-based") is not flux-preserving: emission which cannot 
be uniquely assigned to a cloud is discarded. 

As in HulO, we adopt the standard CPROPS definitions to 
derive the basic physical properties of GMCs in the MAGMA 
cloud catalogue. The cloud radius is defined as R = 1 .91<tr pc, 
where an is the geometric mean of the second moments of the 
emission along the cloud's major and minor axes. The veloc- 
ity dispersion a v is the second moment of the emission distri- 
bution along the velocity axis, which for a Gaussian line pro- 
file is related to the FWHM linewidth, Av, by Av = V81n2cr v . 
The CO luminosity of the cloud Leo is the integrated flux 
scaled by the square of the distance, i.e. 



Leo [K km s" 1 pc 2 ] = D (ST}) SvSxSy 



(1) 



where 7] is the brightness temperature of an individual voxel, 
D is the distance to the LMC in parsecs (taken to be 5 x 10 4 ), 
5x and 5y are the angular pixel dimensions in radians, and Sv 
is the width of one channel in km s -1 . The mass of molec- 
ular gas estimated from the GMCs CO luminosity, M\ um is 
calculated as 



M lum [M Q ]=4.4 



Xco 



2x lO^cm-^Kkms" 1 )- 



rLco , (2) 



where Xco is the assumed CO-to-H2 conversion factor, and 
a factor of 1.36 is included to account for the mass contri- 
bution of helium. Note that the fiducial value of Xco used 
by CPROPS is X C o = 2.0 x 10 20 cnr 2 (Kkms" 1 )" 1 . Al- 
though larger values have been derived for Xco based on 
the assumption of virial equilibrium (Fu08, HulO), in this 
work we do not assume virialization, and so for simplicity 
we calculate M\ am based on this fiducial estimate of Xco- 
iLerov et ail d201 11) use a dust-based method to derive Xco = 
3.0 x 10 20 cm" 2 (Kkms" 1 )" 1 for the LMC, based on com- 
parison of Spitzer IR, NANTEN CO, and ATCA+Parkes H I 
maps. However, they caution that their estimate may be bi- 
ased high by their assumption that the ratio of 160-/im dust 
opacity to gas column density (t^o/A^h) is the same for H I- 
and H2 -dominated regions. There are a number of indications 



8 




Figure 5. Zoom-in of the southeastern portion of the LMC, showing the integrated CO intensity from Figure[T]in grayscale and the sizes and orientations of 
the CO "islands" as black ellipses. The major and minor axes of the ellipses as drawn have been extrapolated to the limit of perfect sensitivity but have not been 
corrected for the finite resolution of the telescope beam. 




EA (degrees; J3000) EA (degrees; J3000) 



Figure 6. Same as Figure|5] but for the (left) northern and (right) western portions of the LMC. 
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Table 4 

Cloud Properties for Islands Decomposition 



ID 


a (2000) 


<5 (2000) 


VTsr 


Maj. 


Min. 


P.A. 


R 




Av 


(5(Av) 


Lco b 


SL 


M vLr c 


SM 


"vox 


T p k 


Isol? 


NT 1 


Hen c 




(h:m:s) 


(":':") 


(km/s) 


(pc) 


(pc) 


(°) 


(pc) 




(km/s) 














(K) 




ID 


ID 


Al 


04:46:50.2 


-67:18:22 


252.6 


10.6 


4.8 


109 


6.3 


0.53 


1.20 


0.41 


8.3 


0.44 


9.4 


1.13 


76 


1.03 


Y 






A2 


04:46:53.0 


-67:13:34 


249.4 


11.2 


7.8 


163 


15.2 


0.42 


1.33 


0.22 


27.6 


0.16 


28.1 


0.65 


273 


1.40 


N 


2 




A3 


04:47:09.6 


-67:12:00 


249.3 


24.4 


5.6 


140 


16.4 


0.23 


1.94 


0.21 


30.2 


0.21 


64.1 


0.53 


295 


1.30 


N 


2 




A4 


04:47:09.8 


-67:09:25 


265.4 


5.9 


5.3 


126 


5.9 


1.33 


0.47 


1.18 


3.6 


1.14 


1.4 


2.68 


32 


0.74 


Y 






A5 


04:47:11.5 


-67:06:58 


261.9 


14.5 


5.8 


97 


13.2 


0.43 


0.59 


0.34 


10.5 


0.26 


4.8 


0.92 


79 


1.44 


Y 


2 




A6 


04:47:26.4 


-69:10:19 


239.8 


28.0 


9.3 


127 


28.6 


0.09 


1.06 


0.10 


80.5 


0.03 


33.4 


0.23 


842 


2.78 


Y 


1 




A7 


04:47:30.2 


-67:13:01 


257.3 


21.0 


9.8 


171 


25.5 


0.20 


2.61 


0.08 


125.5 


0.06 


180.9 


0.28 


1247 


1.91 


N 


2 




A8 


04:47:53.0 


-67:12:50 


265.0 


15.7 


5.3 


146 


12.0 


0.26 


1.08 


0.15 


34.5 


0.10 


14.5 


0.39 


325 


2.09 


N 


2 




A9 


04:47:55.7 


-67:21:04 


256.9 


14.3 


8.1 


20 


18.1 


0.52 


1.72 


0.68 


26.8 


0.45 


55.6 


1.73 


121 


1.22 


N 


2 




A10 


04:48:07.2 


-67:19:30 


258.3 


7.3 


5.8 


1 16 


8.5 


0.87 


1.14 


0.52 


11.3 


0.45 


11.5 


1.45 


68 


1.50 


N 


2 




Al 1 


04:48:10.8 


-67:17:28 


257.8 


5.4 


2.4 


99 






1.00 


0.75 


6.3 


0.49 






41 


1.24 


Y 


2 




A12 


04:48:56.9 


-68:22:23 


250.9 


15.9 


8.9 


162 


20.6 


0.26 


2.11 


0.13 


58.7 


0.16 


95.4 


0.39 


401 


2.13 


N 


5 


N76 


A13 


04:48:57.1 


-69:10:16 


239.9 


23.8 


11.1 


171 


29.4 


0.11 


1.68 


0.08 


130.9 


0.03 


86.1 


0.18 


1249 


2.74 


N 


4 


N77D 



Note. — The full table is available as a machine-readable file. 
"8x is used to denote the fractional uncertainty in x. 
b Units are 10 2 Kkm s _1 pc 2 . 
c Units are 10 3 M Q . 

d Overlapping cloud(s) in Fukui et al. 12008). 
c Associated H II region from Henize ( 1956). 



Table 5 

Cloud Properties for Physical Decomposition 



ID 


a (2000) 


5 (2000) 


W 


Maj. 


Min. 


P.A. 


R 


5R* 


Av 


5(Av) 




SL 


M c 


SM 


"vox 


Tpk 


Isol? 


Isl. 




(h:m:s) 


(":':") 


(km/s) 


(pc) 


(pc) 


(°) 


(pc) 




(km/s) 














(K) 




ID 


Bl 


04:46:50.2 


-67:18:22 


252.6 


10.6 


4.8 


109 


6.3 


0.60 


1.20 


0.38 


8.3 


0.49 


9.4 


1.06 


76 


1.03 


Y 


Al 


B2 


04:46:53.0 


-67:13:34 


249.4 


11.2 


7.8 


163 


15.2 


0.42 


1.33 


0.22 


27.6 


0.16 


28.1 


0.66 


273 


1.40 


N 


A2 


B3 


04:47:09.6 


-67:12:00 


249.3 


24.4 


5.6 


140 


16.4 


0.24 


1.94 


0.19 


30.2 


0.21 


64.1 


0.53 


295 


1.30 


N 


A3 


B4 


04:47:09.8 


-67:09:25 


265.4 


5.9 


5.3 


126 


5.9 


1.34 


0.47 


1.12 


3.6 


1.09 


1.4 


2.69 


32 


0.74 


Y 


A4 


B5 


04:47:11.5 


-67:06:58 


261.9 


14.5 


5.8 


97 


13.2 


0.44 


0.59 


0.33 


10.5 


0.27 


4.8 


0.96 


79 


1.44 


Y 


A5 


B6 


04:47:26.4 


-69:10:19 


239.8 


28.0 


9.3 


127 


28.6 


0.09 


1.06 


0.10 


80.5 


0.04 


33.4 


0.22 


842 


2.78 


Y 


A6 


B7 


04:47:26.6 


-67:11:42 


259.0 


9.3 


4.4 


172 






1.55 


0.28 


37.4 


0.32 






141 


1.91 


N 


A7 


B8 


04:47:30.7 


-67:13:34 


256.3 


14.7 


11.0 


18 


22.6 


0.41 


1.45 


0.26 


75.6 


0.31 


49.1 


0.68 


293 


1.85 


N 


A7 


B9 


04:47:53.0 


-67:12:50 


265.0 


15.7 


5.3 


146 


12.0 


0.25 


1.08 


0.15 


34.5 


0.10 


14.5 


0.39 


325 


2.09 


N 


A8 


B10 


04:47:55.7 


-67:21:04 


256.9 


14.3 


8.1 


20 


18.1 


0.53 


1.72 


0.72 


26.8 


0.44 


55.6 


1.72 


121 


1.22 


N 


A9 


Bll 


04:48:07.2 


-67:19:30 


258.3 


7.3 


5.8 


116 


8.5 


0.86 


1.14 


0.51 


11.3 


0.44 


11.5 


1.52 


68 


1.50 


N 


A10 


B12 


04:48:10.8 


-67:17:28 


257.8 


5.4 


2.4 


99 






1.00 


0.71 


6.3 


0.48 






41 


1.24 


Y 


All 


B13 


04:48:53.8 


-69:10:05 


239.7 


16.5 


7.8 


157 


19.1 


0.25 


1.69 


0.13 


86.4 


0.09 


57.0 


0.40 


431 


2.74 


N 


A13 



Table 6 

Cloud Properties for Data-Based Decomposition 



ID 


a (2000) 


5 (2000) 


v, sr 


Maj. 


Min. 


P.A. 


R 


SR" 


Av 


<5(Av) 


Lco h 


SL 


M vir c 


SM 


"vox 




Isol? 


Isl. 


Phy. 




(h:m:s) 


(»:':") 


(km/s) 


(pc) 


(pc) 


(°) 


(pc) 




(km/s) 














(K) 




ID 


ID 


CI 


04:46:50.2 


-67:18:22 


252.6 


10.6 


4.8 


109 


6.3 


0.61 


1.20 


0.40 


8.3 


0.47 


9.4 


1.10 


76 


1.03 


Y 


Al 


Bl 


C2 


04:46:53.0 


-67:13:34 


249.4 


11.2 


7.8 


163 


15.2 


0.38 


1 .33 


0.21 


27.6 


0.15 


28.1 


0.63 


273 


1.40 


N 


A2 


B2 


C3 


04:47:01.7 


-67:10:52 


247.9 


10.8 


6.7 


136 


13.1 


0.57 


1.20 


0.32 


16.0 


0.45 


19.5 


0.86 


95 


1.30 


N 


A3 


B3 


C4 


04:47:09.8 


-67:09:25 


265.4 


5.9 


5.3 


126 


5.9 


1.39 


0.47 


1.16 


3.6 


1.14 


1.4 


2.69 


32 


0.74 


Y 


A4 


B4 


C5 


04:47:11.5 


-67:06:58 


261.9 


14.5 


5.8 


97 


13.2 


0.40 


0.59 


0.36 


10.5 


0.25 


4.8 


0.95 


79 


1.44 


Y 


A5 


B5 


C6 


04:47:16.6 


-67:13:05 


250.4 


14.0 


5.8 


145 


13.0 


0.58 


1.73 


0.33 


21.0 


0.61 


40.6 


0.90 


100 


1 .05 


N 


A3 


B3 


C7 


04:47:22.8 


-69:09:54 


240.0 


17.8 


7.5 


1 12 


19.2 


0.14 


0.80 


0.13 


52.9 


0.07 


12.9 


0.29 


410 


2.78 


N 


A6 


B6 


C8 


04:47:26.2 


-67:14:13 


255.9 


14.2 


6.9 


169 


16.0 


0.64 


1.66 


0.47 


37.1 


0.82 


46.0 


1.07 


112 


1.45 


N 


A7 


B8 


C9 


04:47:26.6 


-67:11:42 


259.0 


9.3 


4.4 


172 






1.55 


0.27 


37.4 


0.31 






141 


1.91 


N 


A7 


B7 


C10 


04:47:37.4 


-67:12:32 


256.9 


7.6 


5.7 


42 


8.6 


1.17 


1.19 


0.46 


28.3 


0.51 


12.6 


1.58 


80 


1.85 


N 


A7 


B8 


CI 1 


04:47:46.8 


-69:12:25 


238.6 


12.7 


6.0 


9 


12.9 


0.40 


1.07 


0.39 


17.6 


0.44 


15.4 


0.94 


83 


1.77 


N 


A6 


B6 


C12 


04:47:53.0 


-67:12:50 


265.0 


15.7 


5.3 


146 


12.0 


0.24 


1.08 


0.15 


34.5 


0.10 


14.5 


0.41 


325 


2.09 


N 


A8 


B9 


C13 


04:47:55.7 


-67:21:04 


256.9 


14.3 


8.1 


20 


18.1 


0.52 


1.72 


0.71 


26.8 


0.44 


55.6 


1.73 


121 


1.22 


N 


A9 


B10 
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that the FIR emissivity of dust grains is enhanced in regions 
of higher gas density (e.g . jBernard et al.lll999tlStepnik et alj 
I200a iParadis et alj 12009: Pla nck Collaboration et al.1 1201 ll) . 
which would cause the dust-based approach to overestimate 
the amount of molecular gas. 
The virial mass is estimated as 

M vn - [M Q ] = 1040a 2 fl , (3) 

which assumes that molecular clouds are spherical with trun- 
cated p oc r" 1 density profiles (Ma cLaren et al.| [l988). Natu- 
rally, no estimate of M vn - is possible for clouds which are not 
spatially resolved in both dimensions. CPROPS estimates 
the error associated with a cloud property measurement using 
a bootstrapping method, which is described in Section 2.5 of 
RL06. 

As emphasized by RL06, the resolution and sensitivity of 
a dataset influence the derived cloud properties. In order to 
reduce these observational biases, they recommend extrapo- 
lating the cloud property measurements to values that would 
be expected in the limiting case of perfect sensitivity (i.e. a 
brightness temperature threshold of OK), and correcting for 
finite resolution in the spatial and spectral domains by decon- 
volving the telescope beam and width of a spectral channel 
from the measured cloud size and linewidth respectively. The 
procedures that CPROPS uses to apply these corrections are 
described in RL06 (see especially their Figure 2 and Sections 
2.2 and 2.3); briefly stated, a linear extrapolation is used for 
<tr and <j v while a quadratic extrapolation is used for Leo- Un- 
less otherwise noted, the cloud property measurements used 
in the present work have been corrected for resolution and 
sensitivity bias. 

Figures and |6]provide more detailed views of the CO in- 
tensity image with the positions and sizes of the CPROPS "is- 
lands" overlaid. The major and minor axes of the ellipses as 
drawn have been extrapolated to the limit of perfect sensitivity 
but have not been corrected for finite spatial resolution. 

3.3. Catalog Reliability and Completeness 

The overall reliability of the cloud catalogs was checked by 
examining spectra in the CO data cube; in nearly all cases 
there was little doubt that CO emi ssio n had been detected. 
Similarly, as discussed in Section 12.41 we found relatively 
few discrepancies with the NANTEN cloud catalog. How- 
ever, clouds identified as being distinct may still be connected 
at a lower brightness level. We attempted to distinguish "iso- 
lated" islands of emission by excluding islands which connect 
to neighboring islands when contoured down to the 1 .5er level. 
To identify such isolated islands more strictly, we disqualified 
neighboring islands that can be connected through a vertex, 
in addition to islands that connect via a shared face, which is 
the more conventional definition for connected emission. The 
number of isolated islands is tabulated in Table |3j isolated 
clouds are defined in the same way for the other two decom- 
positions. Of course, even the "isolated" regions of emission 
may join with other regions at lower contour levels that we 
are insensitive to. 

Clouds may also be difficult to distinguish because they lie 
along the same line of sight. Because of the low filling frac- 
tion of CO clouds and the face-on aspect of the LMC, we 
believe that chance superpositions of clouds are unlikely. In 
particular, we find very few multi-peaked CO profiles, which 
would be indicative of overlapping clouds at distinct veloci- 
ties. Such an overlap occurs in only 1.6% of positions for the 



"islands" catalog, 2.5% for the "physical" decomposition and 
4.2% for the "data-based" decomposition. Since only strong, 
well-separated components would be identified as distinct, we 
caution that these statistics depend on the signal-to-noise in 
the data cube. 

While the MAGMA maps reveal many previously unknown 
low-luminosity clouds, the exclusion of NANTEN clouds 
with fluxes less than 1 .2 x 10 5 K km s" 1 arcsec 2 from the sam- 
ple definition (corresponding to a luminous mass of M m ; n = 
3 x 10 4 Mq, assuming a Galactic conversion factor) provides 
a minimum completeness limit for the sample. This limit 
corresponds to about twice the completeness limit quoted 
by Fu08 [corrected for a different assumed Xco-factor of 
7 x 10 20 cm~ 2 (Kkm s _1 ) ], suggesting that all clouds with 
mass above M m \ n should have been detected. However, it 
should be kept in mind that both the NANTEN and Mopra 
surveys are brightness temperature limited, so it is not pos- 
sible to determine completeness limits for either survey with 
absolute certainty: clouds covering larger areas and/or having 
larger linewidths may still escape detection. The minimum 
detectable cloud flux (assuming a 3cr brightness temperature 
in at least 3 spectral channels) depends on the cloud size, but 
is roughly 10 2 5 K km s" 1 pc~ 2 for a radius of ~10 pc (com- 
parable to our map resolution). As this corresponds to 1400 
M Q , a factor of 20 below our assumed completeness limit, 
we consider it unlikely that the completeness limit is strongly 
affected by our map sensitivity. We note that while the NAN- 
TEN clouds that form the parent sample of our survey lie 
above the completeness limit, their resolved sub-components 
are likely responsible for populating the regime between the 
detection and completeness limits. 

3.4. Tables of Cloud Properties 

The cloud properties derived by CPROPS are summarized 
in Tables @] [5] and [6] for the "islands" catalog, "physical" 
decomposition, and "data-based" decomposition repectively. 
Only the first 13 lines of each table are shown; full tables are 
accessible from the electronic journal. Each cloud is num- 
bered in order of right ascension and given a prefix of A, 
B, or C to distinguish the three cloud sets. We tabulate in 
columns (l)-(4) the ID number, center of each cloud in de- 
grees of R. A. and Dec. and in LSR velocity (corrected for our 
subtraction of a galactic rotation model). Columns (5)-(7) 
give the major and minor axes of the cloud, corrected for sen- 
sitivity but not resolution effects, and the position angle of the 
major axis (measured in the usual sense from north towards 
east). We then tabulate the size R [col. (8)-(9)], linewidth Av 
[col. (10)-(11)], luminosity L co [col. (12)-(13)], and virial 
mass M v ; r [col. (14)— (15)] of each cloud, along with their re- 
spective fractional uncertainties (a x /x). Columns (16) and 

(17) give the number of significant voxels in the cloud and 
the peak brightness temperature among those voxels. Column 

(18) indicates whether a cloud is "isolated" or not. The re- 
maining columns are different for different tables. In Table|4j 
Column (19) lists cross-identifications with the NANTEN 
catalog (Fu08), and Col umn (20) lists cross-identifications 
with the H II nebulae of iHenizel (119561) . The Henize cross- 
identifications are crude, requiring the center of an island to 
be within 5' of the nebula center (or within the nebula's ra- 
dius if larger than 5'), and should be treated with caution. In 
Table |5j Column (19) indicates the ID number of the island 
in Table [4] which contains the cloud. Similarly, in Table [6j 
Columns (19) and (20) indicate the ID numbers of the islands 
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Figure 7. Cumulative luminous mass spectra for clouds defined with the islands (left), physical (middle), and data-based (right) parameter sets (see Sec. 13.21 
for definitions). Blue solid lines are simple power-law fits to points above a completeness limit of 3 X 10 4 Mq (vertical dashed line). Red dot-dashed lines are 
truncated power-law fits over the same range. 
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Figure 8. Cumulative luminous mass spectra for isolated clouds defined with the islands (left), physical (middle), and data-based (right) parameter sets. Blue 
solid lines are simple power-law fits to points above a completeness limit of 3 X 10 4 Mq (vertical dashed line), or 10 4 Mq for the data-based parameter set. Red 
dot-dashed lines lines are truncated power-law fits. 
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Figure 9. Cumulative virial mass spectra for clouds defined with the islands (left), physical (middle), and data-based (right) parameter sets. Blue solid lines are 
simple power-law fits to points above a completeness limit of 3 X 10 4 Mq (vertical dashed line). Red dot-dashed lines lines are truncated power-law fits. 
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and the "physical" clouds in Tables[4]and[5]which contain the 
cloud. 

3.5. Distribution of Masses 

The differential cloud mass distribution, expressed as the 
number of clouds AMn a logarithmic mass interval d\nM, is 
commonly fitted with a power law, 



dN 
dlnM 



(4) 



where (3 is the usual power law index for the mass distribu- 
tion {dN /dM oc M@). For distributions with (3 > -2, the total 
molecular mass of a galaxy is dominated by massive GMCs, 
while for distributions with /3 < -2, most of the galaxy's 
molecular mass is located in small clouds, and there must be 
a lower mass limit for the molecular clouds or a turnover in 
the mass distribution at low masses in order for the total mass 
to remain finite. In addition, there may be an upper mass limit 
for clouds, beyond which the power law is truncated. 

The cumulative mass distribution, being the integral of the 
mass distribution, has the same slope as the differential mass 
function expressed in logarithmic space, 



, M 

N(M > M) = — 

. M 



0+1 



(5) 



where Mq represents the maximum mass in the sample, which 
generally scales with the total mass in the sample. However, 
for a distribution which is truncated at Mq there is a sharp 
rolloff in the cumulative mass distribution, 



N(M' >M) = N 



M \ 
Mo) 



1 



(6) 



because N(M' > Mp)=0. In this expr essi on, discussed i n mor e 
detail bv lMcKee & Williams! (I 1997b and lRosolowskvl (120051) . 
No is the number of clouds more massive than 2 I < /(/3+1) M , 
where the distribution begins to turn over (for a meaningful 
truncation to exist one expects No S> 1). 

For this study w e focus on the cum ulative mass spectra 
since, as argued by Rosolowskv (2005), these do not require 
the choice of a bin size to generate a histogra m. We ap- 
ply the "error in variables" method employed bv lRosolowskvl 
(2005) to incorporate the uncertainty in the independent vari- 
able (cloud mass) into the functional fits. The method, im- 
plemented in an IDL program provided by E. Rosolowsky, 
uses an iterative maximum-likelihood approach to jointly es- 
timate the "true" sets of M and N (i.e. without measure- 
ment errors) and the best-fit parameters Mo, No, and f3. Fig- 
ure [7] shows the cumulative luminous mass spectra for the 
clouds derived from the three different decomposition meth- 
ods, while Figure H] shows the same distributi ons considering 
only "isolated" clouds, as defined in Section [33] Truncated 
and non-truncated fits, using the functions given in Equa- 
tions (0 and (O, are shown as red and blue curves respec- 
tively. Fitting is only performed above the completeness limit 
of M m i n = 3 x 10 4 M Q shown as a vertical dashed line, except 
for the isolated data-based clouds (Figure [8}:) for which the 
fit was performed above M m ; n = 10 4 M Q in order to allow for 
a sufficient number of points in the fitting. In most cases the 
power-law and truncated fits yield similar slopes, indicating 
that there is little evidence for truncation. Indeed, none of the 
fits satisfies the A^o 3> 1 criterion for a meaningful truncation 
to exist. 



Figures|9]and[T0]show the corresponding cumulative virial 
mass spectra. The virial mass spectra tend to be less steep 
than the luminous mass spectra and are generally consistent 
with a slope of -2. Again, we find little clear evidence for 
a truncation at the high-mass end. The power-law fit results, 
with and without a high-mass truncation, are summarized in 
Table □ 

3.6. Distribution of Radii 

Figure fTTT a) shows the distribution of cloud radii for the 
three different decomposition sets. The NANTEN GMCs of 
Fu08 are also shown as a magenta histogram. There appears 
to be a preferred scale in each distribution corresponding to 
the observational resolution (vertical dashed lines at ~10 pc 
for MAGMA and —40 pc for NANTEN). A straightforward 
interpretation is that clouds near or below the resolution limit 
are incompletely sampled, while clouds larger than the resolu- 
tion limit tend to be resolved into smaller clouds. Figure fTTT b) 
shows the distributions before sensitivity and resolution cor- 
rections are performed. As expected, these distributions show 
sharp cutoffs below the resolution limit. 

The tendency for CO cloud sizes to be close to the observ- 
ing resolu tion has been n oted in a previous study of LMC 
clouds by llsrael et alJ (120031) . and makes it difficult to recover 
the intrinsic properties of molecular clouds from even auto- 
mated decomposition methods. Undoubtedly this is simply 
a refl ection of the hierarch ical structure of molecular clouds 
(e.g.. lFalgarone et aDll991l) . We therefore place greater em- 
phasis on the results derived for the non-decomposed islands, 
since these represent the largest connected CO structures. Al- 
though identification of these structures is still subject to res- 
olution and sensitivity biases (due to blending in crowded re- 
gions and possible failure to detect weak bridging structures), 
the islands span the largest ranges in size, linewidth, and flux, 
making them useful to study trends in cloud properties. 

3.7. Larson-type Relations 

Correlations among cloud properties are frequently dis- 
cussed as indicators of their surface densities and dynami- 
cal state (iBolatto et al] 120081 HulO). We have plotted the 
correlations between size and linewidth, size and luminosity, 
and luminosity and virial mass in Figures [12] Q~3] and [14] re- 
spectively. Each of the three cloud sets is plotted in a sepa- 
rate panel. Red symbols highlight "isolated" clouds, which 
can be uniquely matched to an island defined out to a 1.5cr 
edge. Note, however, that many "isolated" clouds will break 
into multiple substructures upon further decomposition, so the 
number of isolated clouds decreases in going from panel (a) 
to (c) even as the total number of identified clouds increases 
(Table©. 

Following HulO, we fit a linear regression to each panel us- 
ing the 'BCES' (bivariate, correlated errors with intrinsic scat- 
ter) method of Akritas & Bershady (1996), as implemented 
in a FORTRAN program distributed by M. Bershady. This 
method takes into account the uncertainty in each individual 
measurement, as estimated by CPROPS. We have assumed 
that measurement errors in the two plotted variables are un- 
corrected, although some pairs of parameters (e.g., M\ um and 
M v i r ) are expected to have substantial covariance. We plot in 
each panel of Fig. [T2lFT4l a line representing the BCES bisec- 
tor fit, which bisects the results obtained by choosing either 
the X or Y axis as the independent variable. The gray shaded 
region around the fit provides a rough estimate of the uncer- 
tainty in the slope 7 of the bisector fit, as delimited by two 
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Table 7 

Mass Spectrum Fit Parameters 
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Figure 10. Cumulative virial mass spectra for isolated clouds denned with the islands (left), physical (middle), and data-based (right) parameter sets. Blue solid 
lines are simple power-law fits to points above a completeness limit of 3 x 10 4 Mq (vertical dashed line), or 10 4 Mq for the data-based parameter set. Red 
dot-dashed lines lines are truncated power-law fits. 



lines of slope 7-cr 7 and 7 + <r 7 constrained to pass through 
the center of the scatterplot (as defined by the mean values of 
the logarithms of the two parameters). 

For the "islands" decomposition set, the best-fit relation for 
the size-linewidth relation is given by 



log a v = (-0.72 ± 0.05) + (0.80 ± 0.05) logfl 



(7) 



where a v is measured in km s 1 and R is measured in pc. The 
size-luminosity relation can be fit with 



logLco = (1 .00 ± 0. 10) + (2.23 ± 0.08)logfl 



(8) 



where Leo is measured in K km s 1 pc 2 . Finally, the correla- 
tion between CO and virial mass can be expressed as 



logM^ = (-0.28 ± 0. 14) + ( 1 . 15 ± 0.03) logM m 



(9) 



where M\ um is the CO-derived mass assuming a Galactic value 
for the Xco-factor. Figures [T2TlT4l indicate that the slopes of 
the size-linewidth and size-luminosity relations tend to be- 
come steeper and more scattered as one decomposes the CO 
emission further ("physical" and "data-based" sets). Further 
decomposition primarily affects the size distribution, concen- 
trating it near the observational resolution [see also Fig. [TTT a)l 



and resulting in steeper relationships with linewidth and lu- 
minosity. Restricting the fits to "isolated" clouds also leads to 
a slight steepening of the size-linewidth and size-luminosity 
relations, since the largest clouds again tend to be excluded, 
although the differences in slope are at the 1-1. 5cr level. 

HulO obtained quite similar fits to the MAGMA clouds us- 
ing the BCES bisector, although their size-luminosity relation 
is somewhat less steep (7= 1.88 ±0.08) than derived here. 
The HulO analysis differs from the present one in that an ear- 
lier version of the MAGMA cube with somewhat poorer sen- 
sitivity is used, and slightly different criteria are employed 
to identify true emission. HulO also filter out many of the 
weakest clouds based on signal-to-noise considerations, re- 
sulting in few or no clouds with deconvolved linewidths < 1 
km s" 1 or radii <10 pc, whereas the catalogs presented here 
contain substantial numbers of such clouds, due in part to the 
improved sensitivity. 

3.8. Virial Parameter 
Figure[l5]shows the virial parameter, defined as 
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Figure 11. (a) Distributions of the logarithm of the cloud radii, corrected 
for sensitivity and resolution effects, for the three different decomposition 
sets: islands (black), physical (red), and data-based (blue). The NANTEN 
GMCs of Fu08 are also shown as a magenta histogram. Vertical lines indi- 
cate the resolutions of the MAGMA (black) and NANTEN (magenta) surveys 
respectively, (b) Distributions of radii, uncorrected for sensitivity or resolu- 
tion effects, for the three different decomposition sets. 



Table 8 

CO Properties of Candidate YSOs 



GC09 



W08 



Property 


All 


[8.0]<8 


All 


[8.0]<8 


Number of sources 


1172 


247 


1197 


170 


NANTEN observed 


1147 


246 


1141 


164 


NANTEN detected 


669 


191 


367 


60 


MAGMA observed 


722 


196 


421 


67 


MAGMA detected 


552 


171 


259 


46 


MAGMA non-det, NANTEN det. 


99 


15 


75 
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plotted as a function of luminous mass for the three decompo- 
sition sets, with the gray horizontal line representing a v i r =l- 
Under the assumption that CO luminosity faithfully traces 
cloud mass, the condition a v \ r ~ 1 is usually interpreted as be- 
ing gravitationally bound (cf. Bertoldi & McKee 19921), and 
a v ir 3> 1 as being confined by external pressure (which is not 
included in the calculation of a v ir)- The normalization of a v i r 
depends on the appropriate value for Xco, so a value of Xqo 
is often chosen to bring the median value of a v j r close to 1 
(e.g., HulO). In this paper we do not calibrate Xco in this 
way, but simply examine Figure Q3] for evidence of decreas- 
ing a v i r with increasing cloud m ass, as has been pr eviously 
identified in the outer Galaxy by iHever et al.l (1200 lb . How- 
ever, we find no clear tendency for more luminous clouds to 
be more gravitationally bound (with loga v i r ~ 0), although 
many of the "new" clouds that arise from the decomposition 
of larger clouds have values of log a wu < 0. 

3.9. Comparison with YSO Catalogs 

Infrared fluxes and colors of sources in the LMC de- 
rived from Two Micron All-Sky Survey (2MASS) and 
Spitzer Space Telescope observations have been used by 



Whitn ey et al] d2008h (hereafter W08) and iGruendl & Chul 
(2009) (hereafter GC09) to obtain lists of candidate young 
stellar objects (YSOs). Being less sensitive to dust extinction, 
these observations likely trace star formation at earlier stages 
than optical surveys of H II regions and O & B stars. How- 
ever, the masses and ages of the YSOs are highly uncertain, 
as they rely on comparisons with models that are sensitive 
to the d ust geometry and are p rimarily calibrated on Galac- 
tic data (iRobitaille et alj|2007l) . In addition, YSOs can have 
similar broadband colors to asymptotic giant branch (AGB) 
stars, planetary nebulae, and background star-forming galax- 
ies, which will therefore tend to contaminate the candidate 
lists. W08 and GC09 used different sets of criteria to identify 
YSOs and reject likely contaminants, and thus their candidate 
lists, although of similar length, differ substantially. 

We examined both the W08 list of candidate YSOs (1197 
sources) and the GC09 list of "definite" and "probable" YSOs 
(1172 sources) for spatial coincidence with CO emission. In 
practice, this was done by assuming the YSO position was 
precisely known and examining the overlapping pixel in the 
CO map for significant emission. The results are summa- 
rized in Table |8] We used both the NANTEN and Mopra 
maps for identifying coincident CO emission since the NAN- 
TEN map is more spatially complete whereas the Mopra map 
offers better angular resolution. Figure [16] shows the loca- 
tions of the GC09 YSOs relative to the MAGMA CO clouds, 
with red symbols denoting luminous sources (with magni- 
tudes [8.0]<8.0). It is clear that the GC09 sources show 
a high likelihood of being associated with CO emission — 
especially the luminous sources, a large fraction of which 
(97%) were subsequentl y confirmed to have YSO-like spec- 
tra bv lSeale et al.l d2009). In particular, detection rates in the 
NANTEN survey are 58% for the entire GC09 sample and 
78% for the luminous subsample, compared to 32% for the 
entire W08 sample and 37% for the luminous subsample. On 
the other hand, the GC09 list is likely to be biased against 
YSOs that lie outside of GMCs, since the criteria for identify- 
ing YSOs may favor deeply embedded objects and takes into 
account proximity to NANTEN-detected GMCs (although the 
latter is not generally the determining factor). Thus, while the 
GC09 list may contain a larger fraction of bona fide YSOs, the 
W08 list may contain a more representative sample of YSOs 
across a range of evolutionary stages (but see discussion in 
GC09 regarding possible contaminants in the W08 sample). 

CO detection rates when using the MAGMA image are sig- 
nificantly higher than for the NANTEN CO image, which is 
not surprising since the MAGMA observation region is bi- 
ased towards bright CO emission: the fraction of map pix- 
els with detected CO emission is 18% for MAGMA com- 
pared to 9% for NANTEN. It is noteworthy, nonetheless, that 
very few luminous YSOs are found to coincide with CO at 
the NANTEN resolution while being unassociated with CO at 
the improved MAGMA resolution. Only 8% (15/191) of the 
luminous YSOs meet this criterion, whereas the correspond- 
ing fraction is 15% (99/669) for the GC09 YSOs as a whole. 
Given the coarser resolution of the NANTEN data and its bet- 
ter sensitivity to low-brightness emission, this suggests that 
more luminous YSOs tend to be confined within the bright- 
est CO emission whereas less-luminous YSOs are distributed 
both within and around the brightest CO emission. 

Figures [T7l' a)-(d) show the flux distributions of YSOs in the 
NANTEN and MAGMA observation areas and how the prob- 
ability of CO detection varies with 8.0/im magnitude. We 
began by selecting, for the W08 sample, only those sources 
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Figure 12. Size-linewidth relation for the islands (left), physical (middle), and data-based (right) parameter sets. The solid black line is the fitted relati on using 
the BCES bisector, with the gray regions indicating the uncertainty in slope. The mean relations derived for the inner Galaxy by Solomon et al. 1 1987) and for 
nearby galaxies by Bolatto et al. (2008) are shown as a black dotted line and dashed blue line respectively. Larger (red) symbols represent isolated clouds, with 
the remaining clouds shown as small black symbols. 




Figure 13. Size-luminosity relation for the islands (left), physical (middle), and data-based (right) parameter sets. The solid black line is the fitted relation using 
the BCES bisector, with the gray regions indicating the uncertainty in slope. The mean relations derived for the inner Galaxy by Solomon et al. 1 1987) and for 
nearby galaxies by Bolatto et al. (2008) are shown as a black dotted line and dashed blue line respectively. Larger (red) symbols represent isolated clouds. A 
solid green curve represents the estimated sensitivity limit. 




Figure 14. Correlation between luminous and virial mass for the islands (left), physical (middle), and data-based (right) parameter sets. The solid black line 
is the fitted relation using the BCES bisector, w ith the gray regions indicating the uncertainty in slope. The mean relations derived for the inner Galaxy by 
Solomon et al. 1 1987) and for nearby galaxies by Bolatto et al. (2008) are shown as a black dotted line and dashed blue line respectively. Larger (red) symbols 
represent isolated clouds. Gray diagonal lines indicate Xco=0.4, 4, and 40 X 10 20 cm~ 2 (K km s -1 ) , 
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Figure 15. Virial parameter as a function of luminous mass for the islands (left), physical (middle), and data-based (right) parameter sets. Larger (red) symbols 
represent isolated clouds. The gray horizontal line represents a=l. 
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Figure 16. Comparison of CO emission with young stellar populations in the LMC. The jagged gray contour at the periphery of the map shows the region 
surveyed in CO by NANTEN. Multiple rectangular black contours show the regions mapped in CO with Mopra, based on the NANTEN map. Blue contours 
show integrated CO intensity at a level of 1 K km s _I . Filled and open symbols represent "definite" and "probable" YSOs respectively from Graendl & Chu 
(2009), red circles denoting the high-mass sample and green triangles denoting the intermediate-mass sample. 



which are labeled as "highly probable" YSOs and for which 
[8.0] magnitudes were measured (932 sources); for the GC09 
sample all "definite" and "probable" YSOs were considered. 
We then eliminated sources that fell outside the NANTEN or 
MAGMA observation regions and classified the remainder as 
CO detections or non-detections using the appropriate signal 
detection masks from CPROPS (these correspond to the areas 
covered by CO islands before extrapolation). The most lu- 
minous sources, towards the left of the histograms, show the 
highest rate of CO detection, consistent with extreme youth 
and formation within GMCs. The peaks in the distributions 
near magnitude 9 and falloff at fainter magnitudes suggests 
substantial incompleteness in the YSO catalogs for magni- 
tudes [8.0] > 8.0, which may explain the high fraction of CO 
islands that do not harbor any YSO candidates (62% in the 
case of the MAGMA survey, when compared to the GC09 



catalog). 

Figures [T8l a)-(d) show how the probability that a GMC 
contains a YSO depends on GMC luminosity and virial pa- 
rameter. There is a monotonic increase in the probability of 
containing a YSO at higher CO luminosities. The probabil- 
ity of containing a luminous YSO ([8.0]<8.0 mag) also in- 
creases with M\ um . On the other hand, the probability of host- 
ing a YSO is uncorrelated with the virial parameter a vlI . This 
raises doubts about the ability of a wu to diagnose whether a 
cloud can collapse to form stars. On the other hand, clouds 
hosting detectable YSOs may be already undergoing collapse 
motions and/or feedback effects that could modify their line 
widths (and hence virial parameters). 

4. DISCUSSION 
4. 1 . The GMC Mass Spectrum 
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Figure 17. Correspondence between YSO catalogs and CO emission. Histograms show the distribution of 8/^m fluxes for (a) CO from NANTEN, YSO list 
from GC09; (b) CO from MAGMA, YSO list from GC09; (c) CO from NANTEN, YSO list from W08; (d) CO from MAGMA, YSO list from GC09. Only 
YSOs covered by the specific CO survey are shown, with green and red histograms corresponding to detections and non-detections respectively. Vertical dotted 
lines in panels (a) and (b) indicate the separation between "high-mass" and "intermediate-mass" sources for the GC09 sample. The detection fraction is shown in 
the lower panel, with s/N uncertainties in the bin totals indicated. 



The mass spectrum of GMCs in the inn er Galaxy has 
been estimated to have a slope of j3 = —1.5 dSolomon et al.l 

119871; iRosolows kv 2005); A similar slope has been found 

for M31 (IRosolowskyl 120071: iBlitz et al.l [20071) . whereas a 
much steeper slope ((3 = -2.6) has been determined for M33 
(Engargiol a et al] 120031) . In the LMC we find slopes rang- 
ing from (3 = -2.3 to -2.9 for the luminous mass spectrum 
and j3 = -2.0 to -2.6 for the virial mass spectrum, becoming 
steeper (more negative) as one approaches "data-based" de- 
composition parameters. If a steep power-law slope continues 
below our completeness limit, then a significant or even dom- 
inant fraction of the molecular mass may exist in structures 
that we cannot detect individually. Even the fraction of the 
total luminous mass in clouds found in identified structures 
below the adopted completeness limit is substantial (>20%; 
see Table [3), especially for data-based parameters. On the 



other hand, placing too much CO flux in unidentified clouds 
that fall below the completeness limit would v iolat e the inte- 
grated flux measurements presented in Section l2~4l which are 
a factor of ~2 larger than the extrapolated flux in the identi- 
fied MAGMA clouds (Table©. 

For comparison, inferring mass for the NANTEN GMCs 
from CO luminosity, Fu08 measured a slope of j3 = -1.75 
using a completeness limit of 1.4 x 10 4 M Q (corrected to 
our assumed Xco-factor). As noted by Fu08, the fitted mass 
spectrum is quite sensitive to the adopted completeness limit, 
becoming steeper for higher values of M m m- The fact that 
we find a steeper slope using a higher completeness limit of 
^min = 3 x 10 4 M Q is consistent with this trend; if we reduce 
M m j n to 10 4 M then the slope of the luminous mass spectrum 
for islands flattens to -2.0, more similar to th e Fu08 result. 

In addition, as noted by IRosolowskyl d2005l) . blending ef- 
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log(Luminous mass) [M Q ] log(Luminous mass) [M s ] 




Figure 18. Distribution of islands containing (a) YSOs as a function of CO-derived mass; (b) luminous YSOs as a function of CO-derived cloud mass; (c) YSOs 
as a function of virial parameter; (d) luminous YSOs as a function of virial parameter. The upper panels show the distribution of values while the lower panels 
show the fraction of islands containing YSOs, with uncertainties in the bin totals indicated. 



fects will tend to flatten the observed mass spectrum, as sev- 
eral less massive clouds can appear to be a single massive 
cloud. The degree of blending depends on both angular reso- 
lution (which determines whether two objects to be separated) 
and sensitivity (which determines whether two objects can 
be merged at a lower contour level), as well as being sensi- 
tive to the decomposition technique employed. The MAGMA 
maps have better resolution but somewhat poorer brightness 
sensitivity than the NANTEN maps. Both effects will tend 
to separate objects which appear blended in the NANTEN 
maps, steepening the mass spectrum. Of course, even the 
MAGMA data likely suffer from some residual blending: Fig- 
ure |HJ a) suggests a steeper power-law slope (/3 w -2.6 rather 
than -2.3) when islands that connect with other islands at the 
1.5<t contour level are excluded from the fit. 

Even stronger variations in the mass spectrum result when 
islands are decomposed into substructures. Not performing 
decomposition always yields the flattest spectra, while more 
aggressive decomposition tends to depopulate the high-mass 



end and steepen the spectra. There is little doubt that the size 
spectrum following decomposition bears the strong imprint 
of the angular resolution of the observations (Figure fTTI). and 
thus the mass spectrum would be expected to as well. Un- 
less decomposition parameters can be chosen that are tuned 
to size scales that are both physically meaningful (e.g., related 
to a Jeans mass) and well-resolved, it seems doubtful that any 
decomposition method can produce mass spectra that are free 
from observational biases, and the spectrum of "islands" may 
thus prove to be the more useful diagnostic. 

We therefore concu r with previous authors (e.g., She th et aT] 
l2008tlReid et"aT]l2010h that comparisons of mass spectra ob- 
tained with different techniques and from different observa- 
tional data sets must be made with caution. The number and 
properties of the structures identified depend on the dynamic 
range (both spatial and intensity) of the observations, and de- 
ciding which substructures should be considered distinct in- 
troduces additional ambiguity. Choosing an appropriate com- 
pleteness limit is also critically important — too low a com- 
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pleteness limit causes low-mass clouds to be missed, while 
too high a limit leads to sampling and truncation effects at the 
high-mass end. Yet, in the presence of blending, any claim to 
have constructed a complete sample of clouds above a certain 
mass is open to question. Fortunately, the low incidence of 
cloud overlap suggests that line-of-sight blending is not a se- 
vere issue for molecular gas in the LMC, so future maps with 
improved resolution and sensitivity can be expected to mini- 
mize blending. In addition, the use of hierarchical decompo- 
sition methods (e.g. Rosolowsk v et al.120 08; Kauffman n et alj 
120101) will help to characterize the observed structures more 
robustly. 

We note that the distribution of column densities, as shown 
in Figure [4] can be derived without resorting to a cloud de- 
composition, although it is still very dependent on angular res- 
olution: at a scale of 1 1 pc, we are likely smoothing out small- 
scale fluctuations. A lognormal distribution of volume densi- 
ties is expected for isothermal turbulence with a given Mach 
number M, as a result of successive density jumps p\ /po oc 
M 2 produced by shocks (Passot & Vazquez-Semadeni 1998). 
This may translate into a lognormal distribution of column 
densities if the the correlation length for the density field 
is comparable to the cloud si ze, so the largest fluctuations 
occur on the largest scales (Vazquez-Semad eni & Garcfa 
l200l . Detailed mapping of Galactic clouds, however, has 
revealed evidence for departures from the lognormal dis- 
tribution in star-forming clouds, w ith a power-law tail ap- 
pearing at large column densities ( Kainul ainen et al.1 l2009t 
Froebrich & Rowles 2010). In contrast, non-star-forming 
clouds appear to lack these power-law tails. We see no power- 
law excess in our CO-based column density distribution, sug- 
gesting that the scales probed by our observations are still 
dominated by turbulence. In fact, the column density appears 
to be truncated at the high end — a possible indication of opac- 
ity effects in the CO line. We plan to investigate this issue 
in further detail using the 13 CO MAGMA data (Ott et al., in 
preparation). 

4.2. The Larson Scaling Relations 

Correlations between size, density, and line widt h for 
molec ular clouds were discussed in a seminal paper bv fLarsonl 
dl981l) . and have been come to be known as "Larson's laws." 
One of these laws, a roughly linear relationship between 
density and size, implies a constant column density, and 
may be an artifact of the limited range of column densities 
that CO observations are sensitive to, coupled with the ten- 
dency of most lines of sight to sample gas at co lumn densi- 
ties near the observational sensitivity limit ("e.g. JScaTolll99Ct : 
lBallesteros-Paredes~& Mac Low 2002). Consistent with these 
findings, we find the slope of the Lqq-R relation to lie fairly 
close to our estimated sensitivity limit [Figure [T3j a )L with a 
slope for the islands decomposition set that is only slightly 
steeper than the constant column density slope of 2. Fur- 
ther decomposition tends to remove the large R clouds and 
increase the scatter in CO surface brightness, leading to a 
steeper fitted relation [Figure [T3l b)-(c)l. 

On the othe r hand, the size-linewidth relation identified by 
iLarsonl (119811) is not subject to obvious observational biases 
and has frequently been discussed as a real property of in- 
terstellar t urbulence. A relation of the form <r v oc R° 5 was 
derived by Sol omon et al.l dl987l) for the inner Galaxy, and 
a similar rela t ion w ith slope of 0.60 ±0.10 was derived by 
iBolatto eta l. (2008) for extragalactic clouds. Both studies 



have found the normalization of the size-linewidth relation, in 
combination with the mass surface density estimated from CO 
observations, to be consistent with gravitational equilibrium 
(M m 5a 2 R/G). This has been taken as strong evidence that 
molecular clouds are gravitationally bound, although this situ- 
ation does not app ly to lower mass clouds in the outer Galaxy 
dHeveretalJl200l . 

For the MAGMA islands we fit a size-linewidth relation of 
a v oc R &Q±0 05 (Figure [T21 . significantly steeper than previ- 
ous studies. Nonetheless there is considerable scatter around 
the mean relationship: in fact many of the M AGMA is- 
lands are consistent with the Bolatt oet alJ ((2008) relation, al- 
though the majority have smaller linewidths than would be 
predicted for their size. Decomposition of the islands leads 
to a further steepening of the relationship, as large clouds 
are preferentially eliminated from the right side of the di- 
agram. In addition, we find a greatly increased scatter in 
the size-linewidth relation: the decomposed substructures of- 
ten scatter well away from the mean relationship. At face 
value this seems inconsistent with the relation having its ori- 
gin in a turbulent cascade spanning a wide range of scales 
dMac Low & Klessenll2004l) . On the other hand, a CO-based 
decomposition may fail to properly separate structures in 
three-dimensional space that appear superposed in velocity, 
creating additional scatter in the relat ion (see discussion in 
iBallesteros-Paredes & Mac LowH2002l) . While distinct veloc- 
ity components along the line of sight are rarely seen in our 
data, blended velocity profiles would be much more difficult 
to distinguish. 

Figure [14] shows a nearly linear relationship between CO 
luminos ity and virial mass, as has been noted in previous 
studies dSolomon et al.lll987t IBolatto et ai1l2008l) . Although 
our virial masses tend to be higher than our luminous CO 
masses, this may be due to our choice of a Galactic value for 
X C o . Doubling the value of X C o to 4 x 1 20 cm" 2 (K km s" 1 )"' 
would place the two mass scales in reasonable agreeement. 
We are hesitant to use the virial mass to derive a value for 
^co. given that M v ; r oc a 2 R and M\ am oc a v R 2 Tco, so for a 
modest range in CO temperature there is an intrinsic corre- 
lation between M v ; r and M\ um . The scatter in the M vw -M\ um 
relationship, also apparent as a variation in the virial parame- 
ter (Figure[T5l). is quite large, and appears to correlate with the 
scatter in the size-linewidth relation, consistent with the M v - U — 
M\ nm relation ar ising algebraically, as originally discussed by 
Malonev ( 1990). Moreover, Figure Q3] does not reveal a trend 
that more mass ive clouds are more strongly bound, as had 
been found bv lHever et al.l (120011). 

A recent study by iHever et ail d2009l) noted a correla- 
tion between a v /R i) 5 and mass surface density E in the 
Boston University-FCRAO Galactic Ring Survey, with a 
slope that appears consistent with that expected for gravita- 
tionally bound clouds, a v /R 5 = (7rG/5) 1/2 S 1/2 . While such 
a correlation arises t rivially when using virial masses to de- 
rive E, the fact that [Hever et ail d2009l) use masses from an 
LTE analysis of 13 CO leads them to argue in favor of clouds 
being gravitationally bound. However, not unlike the correla- 
tion between M v - n - and M\ um found in this paper, a correlation 
between er,,//? 0,5 and £ arises naturally because both quan- 
tities involve a v /R. The use of dust emission or extinction 
to estimate GMC surface densities may circumvent this is- 
sue. While difficult within the Galactic plane due to line-of- 
sight projection, dust-based approaches applied to the LMC 
will provide an important test of whether the slope and nor- 
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malization of M yn -M\ um correlation are consistent with virial 
equilibrium. 

4.3. Star Formation in LMC GMCs 

It is well known that the correlation between CO emission 
and star formation tracers such as Ha and radio continuum 
emission is rather poor in the LMC (e.g jHughes et al]l2006l) . 
Figure[T]confirms that many luminous H II regions show little 
or no CO emission. On the other hand, Figure [T6l along with 
the detection statistics presented in FiguresTT7land[T8l indicate 
that obscured star formation, as traced by IR-bright YSOs, 
is spatially well-correlated with CO emission. This suggests 
that the differences between Ha and CO distributions can be 
attributed to feedback effects that lead to a rapid dispersal of 
GMCs during the lifetime of the ionizing O & B stars (^10 
Myr). We see possible indications for such evolution in the 
comparison with GC09 YSOs in Figures H3a)-(b): less lu- 
minous (and thus likely longer-lived) YSOs are less likely to 
be associated with CO emission than more luminous YSOs. A 
similar trend is seen for the W08 sources, although with lower 
overall CO detection rates. However, this interpretation is not 
unique, as the less luminous sources may simply lie in clouds 
that are too weak to be detected in CO at our current sensitiv- 
ity. Deeper CO observations are being pursued to investigate 
this possibility 

Recently, Kawamura et al. (2009) compared the locations 
of CO clouds with optical catalogs of young stellar clusters, 
finding that ^2/3 of the <10 Myr old clusters are associated 
with GMCs, the remainder having apparently dispersed their 
natal clouds. This implies a typical duration of ^7 Myr for a 
GMC to host young clusters (although the GMC lifetime may 
be longer if GMCs go through non-cluster-bearing stages). 
Assuming the YSOs detected by Spitzer are much younger 
than 7 Myr, it is thus not surprising that we find nearly all 
of the luminous GC09 YSOs lying within CO clouds. Im- 
portantly, we do not see evidence that any massive GMCs 
lack YSOs, arguing against a long gestation period before 
star formation commences (although the specific rate of star 
formation may vary among GMCs). At present, the limited 
sensitivity of the CO data and the poor characterization of 
YSO ages and masses limit our ability to strongly constrain 
the GMC lifetime, but further work along these lines clearly 
holds tremendous promise for constraining GMC lifetimes, 
thus shedding light on the question of whether GMCs are able 
to achieve gravitational equilibrium. 

5. CONCLUSIONS 

We have presented in this paper a new high-resolution 
imaging survey of CO in the LMC, representing a principal 
outcome of the Magellanic Mopra Assessment (MAGMA). 
The survey targets include the most lu minous GMCs de - 
tected in the second NANTEN survey dFukui et alj 120081) . 
We present three separate catalogs of CO emission struc- 
tures identified using the automated CPROPS algorithm 
dRosolowskv & Lerovll2006h . The "islands" catalog selects 
out all contiguous emission structures, the "physical" catalog 
decomposes the islands into structures designed to resemble 
Galactic GMCs, and the "data-based" catalog decomposes the 
islands into structures to nearly the resolution limit of the sur- 
vey. Based on an analysis of these catalogs, we draw the fol- 
lowing conclusions. 

1. CO "islands" show a steeper mass spectrum than in- 
ferred from previous studies of GMCs in the Galaxy 



dSolomon et al.ll 19871) and of GMCs in the LMC using 
the NANTEN survey. The spectral slope of < -2 sug- 
gests a significant amount of mass in low-luminosity 
clouds. The discrepancy with the earlier NANTEN 
result may reflect uncertainty in the appropriate com- 
pleteness limit above which to fit the spectrum, as well 
as the impact of map resolution and sensitivity in break- 
ing apart blended structures into smaller units. We cau- 
tion that the mass spectrum for decomposed emission is 
very sensitive to the chosen decomposition parameters. 

2. CO "islands" show a roughly constant CO surface 
brightness, as seen in Galactic clouds, but a somewhat 
steeper dependence of line width on radius. Decomposi- 
tion of the islands increases the scatter in both relations 
as large clouds are broken into smaller clouds spanning 
a wide range in linewidth and surface brightness. 

3. The overall distribution of CO column densities in the 
MAGMA field is roughly consistent with a lognormal 
distribution, as expected for isothermal turbulence, al- 
though the distribution may be cut off at the high end 
due to saturation of the CO line. However, these col- 
umn densities are measured over much larger scales (11 
pc) than analogous measurements for Galactic molecu- 
lar clouds. 

4. The nearly linear correlation between virial mass and 
CO luminosity may arise trivially because both quanti- 
ties scale with the product of linewidth and size: virial 
mass scales as a^R and luminosity as a v R 2 Tco- The 
correlation is tightest on scales where the size-linewidth 
relation is tightest, as expected. There is no tendency 
for luminous clouds to show smaller values of the virial 
parameter which would indicate an increasing domi- 
nance of self-gravity. 

5. We find an increased likelihood for more massive 
GMCs to contain YSO candidates, and an increased 
likelihood for more luminous YSO candidates to be as- 
sociated with CO emission. These trends confirm the 
close link between giant molecular clouds (as traced by 
CO emission) and massive star formation: one is rarely 
seen without the other. However, the likelihood of a 
GMC containing a YSO appears unrelated to the virial 
parameter. 

6. We speculate that low-luminosity YSOs may be outliv- 
ing their natal GMCs, based on their lower association 
with CO and tendency to be offset from the brightest 
CO emission. Deeper CO observations and more reli- 
able estimates for YSO masses and ages will be needed 
to confirm this interpretation. 

Future papers in this series will investigate the velocity 
gradients across GMCs, the relationship between atomic and 
molecular gas, and the optical depth of the CO emission (us- 
ing the 13 CO data obtained in parallel). We will also perform 
a comparative analysis of GMCs in the SMC. We expect the 
MAGMA CO maps to have significant legacy value until more 
sensitive and spatially complete mapping surveys are possible 
with multi-pixel receivers, and are releasing the data products 
on web sites hosted by the University of Illinois and CSIRO. 
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